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ABSTRACT 

We  present  high  resolution  (~  0.01  cm'1)  infrared  absorption  spectra  of  the  v4  band  of 
methane  doped  parahydrogen  (CH4/pH2)  solids  produced  by  two  different  techniques:  gas 
condensation  in  an  enclosed  cell  at  T  «  8  K,  and  rapid  vapor  deposition  onto  a  T  «  2  K  substrate 
in  vacuum.  The  spectrum  of  the  rapid  vapor  deposited  solid  contains  a  novel  progression  of 
single  peaks  with  »  5  cm'1  spacing,  superimposed  over  the  known  spectrum  of  CH4  molecules 
trapped  in  sites  of  £>3/ 2  symmetry  in  hexagonal  close-packed  (hep)  solid  pH2.  New  theoretical 
calculations  of  the  rovibrational  transitions  of  a  tetrahedral  molecule  in  an  external  field  of  Oh 
symmetry  permit  the  assignment  of  this  new  progression  to  CH4  molecules  trapped  in  crystalline 
face  centered  cubic  (fee)  regions  of  the  pH2  solid.  Annealing  of  the  rapid  vapor  deposited 
samples  to  T  «  5  K  decreases  the  intensities  of  the  CH4/pH2(fcc)  absorptions,  and  results  in 
intensity  changes  for  parallel  and  perpendicularly  polarized  CH4/pH2(hcp)  transitions.  We 
discuss  these  phenomena,  and  the  narrow  (0.01  to  0.04  cm'1  FWHM)  absorption  linewidths,  in 
terms  of  the  microscopic  structure  of  the  pH,  hosts. 
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I.  INTRODUCTION 

Since  the  discovery  in  the  late  1980’s  by  Oka  and  coworkers1  of  astonishingly  sharp 
infrared  (IR)  absorption  lines  in  doped  cryogenic  parahydrogen  (pH2)  solids,  the  popularity  of 
pH2  as  a  host  for  matrix  isolation  spectroscopy  (MIS)  has  grown  steadily.  The  favorable 
properties  of  pH2  as  a  matrix  host  have  been  discussed  at  length,1"6  and  summarized  quite 
recently,3  so  we  will  not  re-enumerate  them  here.  Rather,  we  focus  our  attention  more  narrowly 
onto  issues  arising  from  the  comparison  of  MIS  experiments  performed  on  pH2  solids  prepared  in 
different  laboratories. 

Liquid  pH2  crystallizes  into  a  hexagonal  close-packed  (hep)  solid  upon  slow  cooling 
below  T  =  13.8  K;  however,  this  process  cannot  generally  be  used  to  prepare  doped  solids.  Two 
principal  methods  for  producing  thick,  transparent,  doped  pH2  solids  have  emerged:  gas 
condensation  in  an  enclosed  cell2'3  held  at  T  «  8  K,  and  rapid  vapor  deposition5’6  onto  a  much 

/■v 

colder  substrate  in  vacuum.  Since  these  solids  can  be  produced  under  highly  nonj-equilibrium 
conditions  they  might  contain  a  variety  of  metastable  microscopic  structures,7  suggesting  a 
corresponding  variety  of  perturbations  to  the  spectroscopy  of  dopant  molecules.  This  manuscript 
addresses  this  connection  between  the  observed  rovibrational  spectroscopy  of  dopant  molecules 
and  the  microscopic  structure  of  the  pH2  host  for  the  specific  case  of  methane  doped 
parahydrogen  (CH4/pH2)  solids. 

IR  absorption  spectra  of  CH4/pH2  solids  produced  by  gas  condensation  in  an  enclosed  cell 
at  T  »  8  K  have  been  previously  reported.8"10  Detailed  analysis  shows  the  existence  of  quantized 
CH4  rotational  states,  with  all  spectral  features  consistent  with  trapping  of  CH4  molecules  in  sites 
of  Dsh  symmetry  (a  variety  of  linear  and  nonlinear  spectroscopies  have  demonstrated  that  pure 
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pH2  solids  grown  in  like  manner  are  polycrystalline  aggregates  of  the  thermodynamically  stable 
hep  phase').  Thus,  these  experimental  and  theoretical  results  provide  a  well  understood 
benchmark  against  which  to  measure  IR  data  obtained  from  other  CH4/pH2  samples. 

Rapid  vapor  deposition  of  pre-cooled  pH2  gas  onto  a  *  2  K  substrate  in  vacuum  can 
produce  millimeters  thick,  optically  transparent  pH2  solids.5  Such  samples  can  be  doped  readily 
by  any  conventional  MIS  dopant  introduction  technique,  and  anecdotal  evidence  suggests  that 
the  lower  condensation  temperatures  result  in  improved  dopant  isolation  efficiency  over  the 
enclosed  cell  method.  However,  the  lower  condensation  temperatures  also  promote  the 
formation  of  metastable  solid  structures,  raising  the  specter  of  multiple  trapping  sites  and 
attendant  inhomogeneous  spectral  perturbations.  IR  and  Raman  spectra  of  pure  pH,  solids 
formed  by  rapid  vapor  deposition  show  that  these  samples:  (a)  are  not  amorphous,  (b)  are  not 
porous,  but  (c)  do  contain  some  non-hep  regions.5  The  Raman  data  indicate  that  the  non-hep 
regions  are  most  likely  facelcenteredTcubic  (fee),5,7  which  differs  from  hep  in  the  stacking 
arrangement  of  close -packed  planes." 

In  what  follows,  we  explore  the  impact  of  these  non-hep  solid  pH2  regions  on  the 
spectroscopy  of  trapped  CH4  molecules.  We  extend  the  theory  of  the  rovibrational  states  of  a 
tetrahedral  molecule  in  a  field  of  Dj/j  symmetry8  (e.g.,CH4  in  hep  pH2)  to  encompass  a 
tetrahedral  molecule  in  an  octahedral  (Ofo)  external  field  (e.gyCH4  in  fee  pH2).  We  show  that  the 
details  of  the  v4  absorption  band  of  rapid  vapor  deposited  CH4/pH2  samples  can  be  explained  as  a 
superposition  of  spectra  for  CH4  trapped  in  both  hep  and  fee  regions. 

This  new  understanding  of  dopant  molecule  spectroscopy  in  return  provides  a  practical 
diagnostic  of  the  microscopic  structures  present  in  solid  pH,  samples  produced  under  variable 
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conditions.  We  hope  the  methodology  developed  during  the  present  collaboration  will  facilitate 
the  comparison  of  MIS  results  obtained  from  doped  pH2  solids  produced  in  other  laboratories. 

II.  THEORY 

A.  fee  Crystal  Field  and  Symmetry  Classification 

We  have  already  developed  the  theory  for  the  analysis  of  the  rovibrational  transitions  of 
methane  in  an  hep  crystal8  based  on  the  extended  group  theory.12  The  same  theoretical 
framework  can  be  applied  to  the  system  of  methane  in  an  fee  crystal,  but  a  modification  is 
necessary  to  take  into  account  the  different  symmetries  between  the  hep  and  fee  crystal 
structures.  In  this  section  we  consider  a  methane  molecule  occupying  a  single  substitutional  site 
of  an  fee  crystal  of  parahydrogen. 

The  structure  of  a  single  substitutional  vacancy  in  an  fee  crystal  belongs  to  a  point  group 
of  Oh,  composed  of  48  symmetry  operations,  while  the  feasible  permutation-inversion  group  of 
methane  is  G24,  which  is  isomorphic  to  a  point  group  of  7j,  composed  of  24  elements.  The 
appropriate  extended  group  Gfcc  for  methane  in  an  fee  crystal  can  be  obtained  by  considering  the 
parity  of  the  whole  system  properly.  The  extended  group  Gfcc  thus  composed  is  isomorphic  to 
the  direct  product  group  of  the  0  and  Trf  point  groups.  The  irreducible  representation  of  the 
group  Gfcc  is  labeled  as  TP  where  T  is  a  representation  of  the  point  group  O,  and  P  is  that  of 
the  group  G24  (T<y).  The  extended  group  for  a  tetrahedral  molecule  in  an  Oh  field  is  given  in 
detail  by  King  and  Homing.13 

The  symmetry  adapted  rotational  wavefunctions  in  terms  of  the  extended  group  Gfcc  are 
obtained  by  considering  the  transformation  of  the  symmetric  top  rotational  wavefunction  |  kJm) 


4 


V 


Final  Draft  27  Apr  99 


under  each  operation  of  the  extended  group.  The  symmetry  adapted  wavefunctions  for  J=  0,  1 , 
and  2  in  the  ground  vibrational  state  are  listed  in  Table  1 .  Here,  the  symmetric  top  rotational 

wavefunction  is  defined  as  | kJm)  =  \lJ  +  l/8/T2  ] 72  ({;£%>})  where  the  phase  of  Wigner’s 

rotation  matrix  D[JJ,  is  the  same  as  before.8,14  We  define  the  crystal  axis  such  that  hydrogen 
molecules  occupy  the  lattice  points  located  at  R,  =  +  /2a2  +  /3a3  where  a,  =  a( i  +  j), 

a2  =  a(  j+  k),  and  a3  =  a(k+  i)  are  the  primitive  vectors  with  the  unit  vectors  i,  j,  k  along  the 
crystal  x,  y,  z- axes.  The  coefficient  a  is  related  to  the  nearest  neighbor  distance  R0  as 
a  =  ,  and  /,  =  0,±1,±2,---  are  the  cell  indices.  To  define  the  molecule-fixed  axis  of 

methane,  one  of  the  C2  symmetry  axes  of  methane  is  taken  as  the  z-axis  and  the  other  two  are 
taken  as  the  x  andy  axes,  so  that  one  of  the  C3  axes  is  oriented  to  point  in  the  (1, 1, 1)  direction  of 
the  molecule-fixed  axis. 

In  order  to  analyze  the  rovibrational  states  of  methane  interacting  with  parahydrogen,  we 
need  to  consider  the  interaction  which  depends  only  on  the  orientation  of  methane.  The 
interaction  potential  for  a  methane  in  an  fee  crystal  has  to  be  invariant  under  any  symmetry 
operation  of  the  extended  group  Gfcc.  Thus,  the  functional  form  of  the  potential  can  be  obtained 
directly  from  the  group  theory.  The  first  anisotropic  potential  is  found  to  be 

V(Q) = *  f  -L  \j%  {DJS  (n)+  d<4>  (n)+  d"\  (n)  +  (n)} 

+  14D«>  fo)  +  5  {d"1  (n)  +  d™  fc>)  +  z>™  (.o)  +  («)}]  (1) 

where  Q  is  the  orientation  of  methane  relative  to  the  z-axis  of  the  crystal.  Note  that  the 
functional  form  of  the  potential  in  Eq.  (1)  depends  on  the  definition  of  the  crystal  and  molecule 
axes.15 
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As  is  seen  in  Eq.  (1),  the  first  anisotropic  interaction  of  methane  in  an  fee  crystal  is 
expressed  in  terms  of  the  fourth  order  rotation  matrix.  On  the  other  hand,  the  lowest  order 
anisotropic  interaction  potential  in  an  hep  crystal  is  found  to  be  expressed  by  the  third  order 
matrix.8  The  difference  comes  from  the  different  symmetry  between  two  crystal  structures:  the 
fee  lattice  is  higher  in  symmetry  than  the  hep  lattice. 

The  difference  is  more  easily  understood  by  constructing  the  interaction  potential  from  all 
the  pairwise  intermolecular  interactions  between  the  methane  and  hydrogen  molecules  in  a 
crystal.  The  first  two  symmetry  adapted  interaction  functions  for  a  single  pair  of  a  methane  and  a 
hydrogen  molecule  are  given  as: 

v3  (R,  to)  =  a 3  {R)-j=  [.Dgj  {co)~  D%  (©)]  (2) 

and 

v4  (R,a>)=  a4  {R)^  [V5  '{d<J  (co)+  D%  (*>)}+  (to)  (3) 

where  R  is  the  distance  between  the  methane  and  hydrogen  molecules  and  co  is  the  orientation  of 
methane  relative  to  the  pair  axis  between  the  two  molecules.  The  symbols  az(R)  and  a4(i?)  are  the 
distance  dependent  coefficients  for  these  interactions.  The  pairwise  potentials  in  Eqs.  (2)  and  (3) 
are  more  concisely  written  as: 

vl{R,co)=al{R)YJA[^ '>(<»)  (4) 

k=-l 

where  1=3,  4  and  A[  is  a  proper  coefficient,  i.e.  Al2  =±i/4 2  ,  A*4  =  V30/l2 ,  and  A\  =  V2l/6 . 
Since  the  interaction  between  methane  and  parahydrogen  is  mainly  due  to  the  dispersion  force 
which  is  nearly  additive,  the  interaction  potential  in  the  crystal  can  be  obtained  as  a  summation 
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of  pairwise  intermolecular  interactions  of  Eqs.  (2)  and  (3).  Using  the  fact  that  the  coefficients 
a3(R)  and  a4(i?)  for  the  dispersion  interaction  are  proportional  to  R'7  and  i?'8,  respectively,16  the  /- 
th  order  interaction  potential  in  the  crystal  V,(Q)  is  expressed  as: 

(5). 

m=-l  k=-l 

Here,  the  symbol  S^m  is  a  crystal  sum  defined  as: 

C=£M0£ C*(0  w 

7=2 

where  <5=7  for  1=3  and  <5=8  for  1=4,  R0  is  the  nearest  neighbor  distance,  Rj  is  the  distance  between 
the  methane  and  the y-th  hydrogen  molecule,  and  the  angle  Q,  .  is  the  orientation  of  the  pair 
axis  between  the  methane  and  the  y-th  hydrogen  molecule  relative  to  the  crystal  axis.  Note  that 
the  crystal  sum  S^m  depends  only  on  the  configuration  of  the  molecules  in  the  lattice. 

In  the  case  of  an  hep  crystal  the  crystal  sum  for  1=3,  m=± 3  is  not  zero,  and  thus  the  third 
order  potential  caused  by  the  pairwise  interaction  of  Eq.  (2)  remains  (see  Eq.  (17)  in  Ref.  8).  On 
the  contrary,  in  the  case  of  an  fee  crystal,  the  first  non- vanishing  crystal  sum  is  that  of  l  =  4 

(m  =  0,  ±4),  namely,  S$  =+A/l4/55s(j4±)4  =-1.5729.  Thus,  the  third  order  potential  does  not  exist 

in  this  case  but  the  fourth  order  potential  resulting  from  the  interaction  of  Eq.  (3)  remains  with 
the  functional  form  given  in  Eq.  (1). 

B.  Order  Estimation 

Since  the  potential  in  an  fee  crystal  is  a  function  of  the  fourth  order  Wigner’s  matrix,  the 
splittings  of  rotational  M  sublevels  are  expected  to  be  small  compared  with  that  in  hep  crystal. 
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The  coefficient  efc  in  Eq.  (1)  is  related  to  the  coefficient  a4  in  Eq.  (3)  as: 

=6S<>4(*.)A/2i 

By  use  of  the  relation15 

3^35  ACch  1 

a^R)  =  -W~I^laAR) 


(7). 


(8) 


we  obtain  the  following  equation  relating  the  coefficient  s{cc  in  Eq.  (1)  and  the  coefficient 

s*ccp  of  the  crystal  field  in  an  hep  crystal  defined  in  our  previous  paper,8 

_,cc  3^30  C/cc  ACCHj  1 
4  “  28/  Ac h4  R0  ^ 

Here,  S?}hcp ,  ACHj  ,  and  ACCH4  are  the  crystal  sums  of  the  hep  lattice,  dipole-quadrupole 
polarizability,  and  quadrupole  polarizability  of  methane,  respectively.16  The  experimental  or 
theoretical  values  are  S(73jhcp  =  0.529/  (where  /  s  V-T ,  [ref.  8]),  R0  =  3.783  A  [ref.  8], 

ACch  /Ach  =  -0.52  A  [refs.  17,  18],  and  s^p  =  -25.8  cm'1  [ref.  8].  Thus,  we  arrive  at  an 
estimated  value  of  s{cc  =  -6.2  cm'1. 

The  rovibrational  states  of  methane  in  an  fee  crystal  are  obtained  by  solving  the  standard 
vibration-rotation  Hamiltonian  matrix  under  the  potential  of  Eq.  (1).  In  figure  1  the  rovibrational 
energy  levels  of  methane  and  the  optical  selection  rules  are  compared  between  the  isolated 
molecule  and  the  molecule  in  the  fee  crystal.  Those  of  the  molecule  in  the  hep  crystal  are  also 
shown  in  figure  1  for  comparison.  The  solid  arrows  in  the  diagram  for  the  isolated  molecule  and 
the  molecule  in  the  fee  crystal  indicate  allowed  transitions.  Since  the  space  fixed  components  of 

the  dipole  moment  belong  to  the  Fl  A2  symmetry  species  in  the  extended  group  Gfcc,  all 
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transitions  become  allowed  for  any  orientation  of  polarization  of  light  relative  to  the  crystal  axis. 
On  the  contrary,  in  the  case  of  the  hep  crystal,  the  parallel  (solid  arrows  in  figure  1)  and 
perpendicular  (dashed  arrows)  transitions  show  different  selection  rules  since  the  hep  crystal  c- 
axis  is  different  from  the  other  a-  and  b-  axes. 

The  J=  1  ground  vibrational  state  belongs  to  Fi  Fx  symmetry  in  the  extended  group. 

Thus,  the  crystal  field  of  Oh  symmetry  does  not  lift  the  degeneracy  of  the  M-sublevels  of  J=l  as 
shown  in  figure  1.  This  is  in  contrast  to  the  case  of  methane  in  the  hep  crystal  where  the  M- 
sublevels  of  J=  1  split  into  two  due  to  the  crystal  field  of  D^h  symmetry.  In  the  same  way,  no 
splitting  of  the  M-sublevels  of  J- 1  is  expected  in  the  excited  vibrational  state  in  the  fee  crystal. 
The  splitting  caused  by  the  crystal  field  of  Oh  symmetry  first  appears  at  v=l,  J=2,  R=  1  levels 

having  F:F2  and  EF2  symmetries.  By  the  use  of  the  crystal  field  value  of  s{cc  =  -6.2  cm'1,  the 
splitting  width  of  this  level  is  numerically  calculated  to  be  ~  0.01  cm'1. 

Thus,  for  the  infrared  absorption  spectrum  of  methane  in  the  fee  crystal  we  expect  only 
one  absorption  line  for  each  rotational  branch  P(l),  Q(l),  and  R(0),  while  the  R(l)  branch  should 
appear  as  a  doublet  with  a  splitting  of  ~  0.01  cm'1.  These  predictions  are  based  on  the 
assumption  that  the  crystal  structure  is  perfect.  If  there  are  some  imperfection  of  the  crystal 
stacking  (as  will  be  discussed  later),the  symmetry  of  the  crystal  field  is  further  reduced,  which 
may  cause  additional  small  splittings  even  in  J=l  rotational  levels. 


III.  EXPERIMENTAL 
A.  Kyoto  University 

Figure  2  shows  a  schematic  of  the  experimental  apparatus  residing  at  Kyoto  University 
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0.018, 0.020,  and  0.024  cm'1,  respectively.  The  slightly  broader  width  for  the  R(l)  peak  may 
reflect  a  contribution  from  the  ~  0.01  cm'1  crystal  field  induced  splitting  predicted  above  at  the 
end  of  section  II.B.  The  observation  of  only  one  absorption  line  for  each  rotational  branch  of 
CH4/pH2(fcc)  indicates  that  the  symmetry  of  the  crystal  field  is  almost  perfectly  Oh- 

This  set  of  four  observed  transitions  is  not  sufficient  to  determine  molecular  parameters 
for  CH4/pH2(fcc).  All  of  these  transitions  satisfy  a  selection  rule  of  AR  =  0  as  is  seen  in  figure  1 . 
However,  without  A R  ^  0  transitions,  the  rotational  constant  and  the  Coriolis  parameter  cannot 
be  determined  independently.  Moreover,  transitions  to  higher  J  levels  are  indispensable  in 
determining  the  crystal  field  parameter  s^cc  quantitatively,  since  the  crystal  field  potential  is  a 
function  of  only  the  fourth  order  Wigner’s  matrix.  Prospective  observations  of  Raman 
transitions,  which  obey  a  A /  =  0,  ±1,  and  +2  selection  rule,  may  allow  us  to  determine  molecular 
parameters  for  CH4/pH2(fcc). 

Transitions  of  CH4/pH2(hcp)  in  figure  5  are  clearly  separated  under  high  resolution,  as 
predicted  by  theory.8,9  Especially  noteworthy  is  a  clearly  resolved  doublet  at  1312.92  cm 1  an 

1312.95  cm'1  which  corresponds  to  the  transitions  (J-2,  R  -1,  M-  0,  AXF2)<-  (J”=l,  M”=0 
A2JX )  and  (J-2,  R  -1,  M’=±  1,  EJX )  <-  (J”=  1,  M”=±  1,  eT2  ),  respectively  (see  Table  III  ol 
Ref.  9).  The  high-resolution  data  in  Table2  allows  us  to  modify  the  molecular  parameters  of 
CH4/pH2(hcp)  in  Ref.  9,  which  will  be  published  separately  along  with  those  of  all  the 
isotopomers  of  methane. 

We  also  observe  small,  but  experimentally  significant,  differences  in  the  transition 
frequencies  between  the  rapid  vapor  deposited  samples  and  the  sample  condensed  in  an  enclosed 
cell,  as  is  seen  in  figure  5  and  Table  II.  These  spectral  differences  are  mainly  explained  by 
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C.  Other  Possible  Microscopic  Structures  for  Vapor  Deposited  pH2  Solids 

Although  we  have  discussed  the  spectra  obtained  from  rapid  vapor  deposited  CH4/pH2 
solids  in  terms  of  trapping  of  CH4  molecules  in  hep  and  fee  crystallites,  it  is  important  to 
consider  whether  other  solid  pH2  structures  might  also  be  consistent  with  our  data.  The  hep  and 
fee  structures  are  only  two  of  the  infinite  number  of  densest  close-packed  structures  derivable  by 
stacking  of  close-packed  planes.11  The  fee  structure  (...ABCABCABC...)  is  actually  the  simplest 
in  that  each  close-packed  plane  has  the  same  translation  relative  to  the  plane  below,  whereas  in 
hep  (...ABABAB...)  each  successive  plane  reverses  its  predecessor’s  lateral  displacement.20 
Patterns  with  longer  repeat  units  (polytypism),  non-periodic  structures  such  as  twin  and  stacking 
faults,  and  even  random  stacked  close-packed  (rep)  structures  are  also  possible,  in  principle.11,20 
However,  if  only  nearest  neighbor  interactions  are  considered,  then  there  are  just  two  types  of 
single  substitutional  trapping  sites  for  any  of  these  structures:  sites  in  layers  (B)  embedded  in 
regions  of  local  hep  stacking  (XABAX,  “hep-like”),  or  local  fee  stacking  (XABCX,  “fcc-like”). 

We  thus  consider  two  extreme  ideal  structures  for  the  transparent  rapid  vapor  deposited 
pH2  solids:  (a)  an  aggregate  of  distinct  perfect  fee  and  hep  crystallites  in  which  ordered  stacking 
sequences  persist  over  dimensions  larger  than  optical  wavelengths,  and  (b)  a  random  stacked 
close-packed  solid  for  which  fcc-like  and  hep-like  structures  are  intermingled  at  the  molecular 
level.  Without  some  independent  knowledge  of  how  CH4  trapping  site  structures  (and  hence 
spectral  perturbations)  should  differ  in  hep  vs.  hep-like  pH2  solids,  we  cannot  make  an 
unequivocal  distinction  based  on  the  present  spectroscopic  data. 

In  an  attempt  to  directly  observe  the  effects  of  increasing  structural  disorder  on  the 
spectroscopy  of  CH4/pH2  solids,  we  performed  a  series  of  sample  in  vacuum  depositions 
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employing  successively  slower  pH2  inlet  flow  rates.  The  slower  flow  rates  result  in  lower 
deposition  temperatures,  which  we  hoped  would  further  enhance  the  formation  of  metastable 
structures.  The  results  are  summarized  in  figure  6,  which  shows  the  Q(l)  region  of  spectra  of 
four  as-deposited  CH4/pH2  samples.  Decreasing  the  pH2  inlet  flow  from  190  to  45  mmol/h, 
traces  (6a)  and  (6b),  indeed  results  in  increased  intensity  for  CH4/pH2(fcc)  absorptions,  in  slightly 
broader  linewidths,  and  in  slightly  improved  alignment  of  the  hep  c-axes  with  the  substrate 
surface  normal.  However,  further  decreasing  the  pH2  inlet  flow  to  1 0  and  then  to  4  mmol/h, 
traces  (6c)  and  (6d),  results  in  relatively  weaker  CH4/pH2(fcc)  absorptions  and  in  narrower  lines; 
the  trend  towards  increased  alignment  between  hep  c-axes  and  the  substrate  surface  normal  is 
continued.  Surprisingly,  all  of  these  samples  exhibited  excellent  optical  transparency,  even  the 
two  deposited  at  the  slowest  flow  rates. 

Absent  a  qualitative  spectroscopic  signature  for  the  formation  of  a  microscopically 
disordered  random  close  packed  solid,  we  conclude  that  the  mixed  hep/fee  crystallite  model  is 
capable  of  explaining  all  of  our  results  to  date.  We  await  the  application  of  theoretical 
techniques  such  as  quantum  Monte  Carlo  simulations  of  trapping  site  structures  in  imperfect 
crystals,  and  the  accumulation  of  a  broader  spectroscopic  database  on  doped  pH2  solids  to 
definititively  resolve  this  issue. 

V.  CONCLUSIONS 

Two  methods  for  preparing  millimeters  thick,  optically  transparent,  doped  pH2  solids 
have  been  presented:  gas  condensation  in  an  enclosed  cell,  and  rapid  vapor  deposition  onto  a 
substrate  in  vacuum.  Rapid  vapor  deposition  of  precooled  pH2  gas  onto  a»2K  substrate  in 
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TABLES 

TABLE  I.  Symmetry  adapted  rotational  wavefunctions  of  the  ground  vibrational  state  expressed 
as  linear  combinations  of  the  symmetric  top  rotational  wavefunction  |  kJm)  and  their 
representation  in  the  extended  group  G. 


J=  0  Ax  Ax 

J=  i  fxfx 

J=  2  EE 


EF2 

f2e 

f2T2 


|0,0,0) 

|0,1,0),  |±1,1j0>  ,  |o,i,±i) ,  I ±1,1-1),  I ±1,1,1) 

\  I  -  2,2  -2) + 1 2,2-2)  + 1  -  2,2,2)  + 1 2,2,2)],  -±=  1 0,2, -2)  + 1 0,2,2)] , 
-j=[-2,2,0)  +  |2,2,0)],  |0,2,0) 

~  | _  2, 2, -2)  - 1 2,2, -2)  + 1  -  2,2,2)  - 1 2,2,2)] ,  -j=[-2,2,0)-|2,2,0)], 

|  ±1,2,0) ,  |  ±1,2,2) 

~ lj ~ 2,2, -2)  + 1 2, 2,-2) - 1 - 2,2,2) - 1 2,2,2)] ,  -j=|o,2,-2)-|0,2,2)], 

-^(-2,2,±1)  +  |2,2,±1)],  |0,2,±l) 

~ [j ~ 2,2, -2) - 1 2, 2, -2) - 1 - 2,2,2)  + 1 2,2,2)] ,  j=j±l,2,-2)-|±l,2,2)], 
-j=[j-2,2,±l)-|2,2,±l)],  |±  1,2,— l) ,  | ±1,2,1) 
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** 


•t 


TABLE  II.  Absorption  peaks  (cm'1),  and  spectroscopic  assignments  for  v4  fundamental  band  of 
CH4/pH2.  The  labels  (5a),  (5b),  and  (5c)  refer  to  the  spectral  traces  and  sample  descriptions  in 
figure  5. 


(5c) _ assignment  (T.  R'.  M'.  Pi  4-  (J''.  M".  F') 


(5a) 

_ (5b) 

1298.236 

1298.240 

1298.312 

1298.892 

1298.898 

1302.943 

1302.943 

1303.331 

1303.332 

1303.598 

1303.600 

1303.846 

1303.849 

1308.259 

1308.255 

1308.377 

130.8.374 

1308.412 

1308.412 

1309.188 

1309.185 

1312.289 

1312.290 

1312.92  sh 

1312.93  sh 

1312.960 

1312.945 

1313.054 

1313.578 

1313.579 

1316.504 

1316.515 

1298.237 

P(l)  hep  || 

P(l)  fee 

1298.894 

P(l)  hep  1 

1302.944 

Q(l)  hep  1 
Q(l)  fee 

1303.600 

Q(l)  hep  || 

1303.848 

Q(l)  hep  1 

1308.261 

R(0)  hep  1 

1308.376 

R(0)  hep  || 
R(0)  fee 

1309.191 

R(0)  hep  J. 

1312.292 

R(l)  hep  1 

1312.923 

R(l)hcp  || 

1312.947 

R(l)  hep  || 
R(l)  fee 
R(l)  fee 

1313.578 

R(l)  hep  1 

1316.510 

R(l)  hep  1 

(0,1,0,  AXF\)  4- (1,0  ,44) 
(0,1, x, 4^)  4-(l ,x,44) 
(0,1,0,44)  4-(l,±l,£4) 

(U,±i,£4)  4- (1,0,44) 

(1,1*44)  <-(U,44) 

(U,±i,E4)  4- (i, ±\,eT2) 

(i,i,o,44)  <-(i,±i,£4) 

( 1 ,0,0,  EAX )  4—  (0,0, 4  4 ) 
(i,o^i,44)4-  (0,0,44) 
(1,0*44)  4-(0,x,44) 
(1,2,±1  ,ET2)  4-(l,±l ,eT2) 

(2,1,±1  ,EJX)  4- (1,0,44) 

(2,1,0,44^)  <-(i,o,4^) 

(2,1,±1,E4)  4-(l,±l,E4) 
(2,1,0, EF\ )  4-(1,x,44) 
(2,i,±i,44)4-(i,x,44) 

(2,1,0,44)  <-(i,±i,£4) 

(2,1, ±2,  £4  )  4-  (1,±1,£4  ) 


||  and  .I  -  polarizations  relative  to  c-axis  in  hep  pH2 
x  -  all  M  substates  remain  degenerate 
sh  -  shoulder 
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FIGURE  CAPTIONS 

FIG.  1.  Rotational  levels  of  the  ground  and  v4  excited  vibrational  states  of:  (1)  isolated  CH4,  (2) 
CH4  in  fee  pH2,  and  (3)  CH4  in  hep  pH2.  Allowed  transitions  are  shown  as  vertical  solid  arrows 
in  the  diagram  of  isolated  andjfcc.  Vertical  solid  and  dashed  arrows  in  the  diagram  of  hep  shows 
the  parallel  and  perpendicular  transitions,  respectively.  The  irreducible  representation  of  each 
rotational  level  in  terms  of  the  extended  group  G  is  listed  at  the  right-hand-side  of  the  energy 
levels  in  each  diagram. 

FIG.  2.  Experimental  schematic  (Kyoto  University). 

FIG.  3.  Experimental  schematic  (Edwards  Air  Force  Base). 

FIG.  4.  Overview  of  the  v4  region  in  IR  absorption  spectra  of  a  rapid  vapor  deposited  CH4/pH2 
sample.  Trace  (a)  is  for  the  as-deposited  sample  at  Tsub  =  2.2  K,  (b)  warmed  to  Tsub  =  4.4  K,  (c) 
cooled  to  Tsub  =  3.0  K,  and  (d)  shows  the  annealed  sample  cooled  to  Tsub  =  2.2  K.  The  sample 
remains  at  each  temperature  for  1  h  as  200  interferometer  scans  (0.1  cm'1  resolution)  are 
accumulated.  During  the  25  minute  long  deposition  the  pH2  inlet  rate  is  190  mmol/h  (sample 
thickness  growth  rate  =  50  pm/min)  and  Tsub  rises  from  2.2  to  3.2  K.  The  final  CH4  concentration 
is  estimated  as  1 10  PPM,  and  the  sample  thickness  as  1.2  mm.  Most  peaks  show  instrument 
limited  linewidths  of  ~  0. 1  cm  .  The  peaks  in  frace  (a)  which  dramatically  decrease  in  intensity 
upon  annealing  are  marked  with  asterisks. 
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FIG.  5  High  resolution  detail  of  the  v4  absorption  region  of  CH4/pH2  samples  produced  by  rapid 
vapor  deposition,  and  by  gas  condensation  in  an  enclosed  cell.  Trace  (a)  shows  the  rapid  vapor 
deposited  sample  at  Tsub  =  2.4  K  immediately  after  deposition.  Trace  (b)  shows  the  rapid  vapor 
deposited  sample  at  Tsub  =  4.5  K,  after  a  double  annealing  cycle:  2.4  K  — >  4.8  K  —>  3.6  K  -»  2.4 
K  -»  4.5  K.  The  sample  remains  at  each  temperature  for  1  h  as  40  interferometer  scans  (0.008 
cm'1  resolution)  are  accumulated.  During  the  1  h  long  rapid  vapor  deposition  the  pH2  inlet  rate  is 
190  mmol/h  (sample  thickness  growth  rate  =  50  pm/min)  and  Tsub  rises  from  2.2  to  3.2  K.  The 
final  CH4  concentration  is  estimated  as  15  PPM,  and  the  sample  thickness  as  2.9  mm.  The  peaks 
in  trace  (a)  which  dramatically  decrease  in  intensity  upon  annealing  are  marked  with  asterisks. 
Trace  (c)  shows  a  spectrum  (0.01  cm'1  resolution)  of  the  enclosed  cell  condensed  sample  cooled 
to  4.5  K  following  the  2  h  long  sample  growth  process  at  T  ~  8  K;  the  isolated  CH4  concentration 
is  estimated  as  2  PPM.  The  symbols  ||  and  _L  indicate  the  predicted  polarizations,  relative  to  the 
hep  crystal  c-axis,  for  allowed  transitions  of  CH4  molecules  trapped  in  hep  pH2. 

FIG.  6  High  resolution  spectra  (0.005  cm'1  resolution)  of  the  Q(l)  feature  for  as-deposited 
CH4/pH2  samples  produced  at  various  pH2  flow  rates.  Details  for  each  1  h  long  deposition:  (a) 
pH2  inlet  rate  is  190  mmol/h,  T^  rises  from  2.2  to  3.2  K,  final  sample  thickness  d  =  2.9  mm,  and 
the  CH4  concentration  CCH4 «  15  PPM,  (b)  pH2  inlet  rate  =  45  mmol/h,  Tsub  rises  from  2.2  to  2.4 
K,  d  =  0.8  mm,  and  CCH4 «  55  PPM,  (c)  pH2  inlet  rate  is  10  mmol/h,  Tsub  rises  from  2.2  to  2.3  K,  d 
-  0.2  mm,  and  CCH4  *  200  PPM,  and  (d)  pH2  inlet  rate  is  4  mmol/h,  Tsub  remains  at  2.2  K,  d  = 

0.08  mm,  and  CCH4  *  550  PPM.  The  symbols  *,  ||,  and  _L  are  defined  in  the  captions  to  figures  4 
and  5. 
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